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ABSTRACT 

We  present  the  early  UV  and  optical  light  curve  of  Type  IIP  supernova  (SN)  2010aq  at  z  =  0.0862,  and  compare 
it  to  analytical  models  for  thermal  emission  following  SN  shock  breakout  in  a  red  supergiant  star.  SN  2010aq  was 
discovered  in  joint  monitoring  between  the  Galaxy  Evolution  Explorer  {GALEX)  Time  Domain  Survey  (TDS)  in 
the  NUV  and  the  Pan-STARRSl  Medium  Deep  Survey  (PS1  MDS)  in  the  g,  r,  i,  and  z.  bands.  The  GALEX  and 
Pan-STARRS  1  observations  detect  the  SN  less  than  1  day  after  the  shock  breakout,  measure  a  diluted  blackbody 
temperature  of  31,000  ±  6000  K  1  day  later,  and  follow  the  rise  in  the  UV/optical  light  curve  over  the  next  2  days 
caused  by  the  expansion  and  cooling  of  the  SN  ejecta.  The  high  signal-to-noise  ratio  of  the  simultaneous  UV  and 
optical  photometry  allows  us  to  fit  for  a  progenitor  star  radius  of  700  ±  200 RQ,  the  size  of  a  red  supergiant  star.  An 
excess  in  UV  emission  two  weeks  after  shock  breakout  compared  with  SNe  well  fitted  by  model  atmosphere-code 
synthetic  spectra  with  solar  metallicity  is  best  explained  by  suppressed  line  blanketing  due  to  a  lower  metallicity 
progenitor  star  in  SN  2010aq.  Continued  monitoring  of  PS1  MDS  fields  by  the  GALEX  TDS  will  increase  the 
sample  of  early  UV  detections  of  Type  II  SNe  by  an  order  of  magnitude  and  probe  the  diversity  of  SN  progenitor 
star  properties. 

Key  words:  supernovae:  individual  (SN  2010aq)  -  surveys  -  ultraviolet:  general 


1.  INTRODUCTION 

Shock  breakout  in  a  core-collapse  supernova  (SN)  marks  the 
first  escape  of  radiation  from  the  blast  wave  that  breaks  through 
the  surface  of  the  star  and  launches  the  SN  ejecta,  heating  it  to 
temperatures  of  105-106  K  (Colgate  1974;  Falk  1978;  Klein  & 
Chevalier  1978;  Matzner  &  McKee  1999).  The  duration  of  this 
radiative  precursor  in  the  UV /X-rays  is  smeared  out  by  light- 
travel  time  effects  to  fit  =  R„/c  ~  2  s  (R*/RG)>  where  R+  is 
the  pre-SN  radius  of  the  progenitor  star.  In  the  presence  of  a 
dense  wind,  or  if  the  shock  breakout  is  aspherical  or  very  weak, 
the  observed  duration  can  be  much  longer  (Campana  et  al.  2006; 
Chevalier  &  Fransson  2008;  Couch  et  al.  2009;  Dessart  et  al. 
2010). 

Following  shock  breakout,  the  early  SN  light  curve  is  powered 
by  radiation  from  the  expanding  and  cooling  SN  ejecta.  In 
hydrogen-rich  Type  II-plateau  SNe  (SNe  IIP),  a  receding  H 
recombination  wave  prolongs  the  optical  brightness  of  the  SN 
until  radioactive  heating  plays  a  dominant  role.  The  luminosity 
and  temperature  evolution  of  the  early  thermal  expansion  phase 
can  be  used  as  a  diagnostic  for  the  radius  of  the  progenitor  star, 
as  well  as  the  explosion  energy  to  ejecta  mass  ratio  (Waxman 
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et  al.  2007).  This  is  an  important  tool,  since  the  direct  detection 
of  SN  progenitors  is  incredibly  difficult,  and  has  only  been 
possible  for  a  small  number  of  nearby  SNe  with  pre-explosion 
high  resolution  imaging  (Smartt  et  al.  2009). 

Recent  observational  advances  have  enabled  the  detection  of 
early  UV  emission  from  core-collapse  SNe  from  hours  to  a  cou¬ 
ple  of  days  after  explosion.  The  Swift  satellite  serendipitously 
detected  luminous  X-ray  bursts  with  X-Ray  Telescope  (XRT) 
from  Type  Ib/c  SNe  2006aj  and  2008D,  which  were  modeled 
by  shock  breakout  in  the  dense  wind  surrounding  a  Wolf-Rayet 
(W-R)  star  (Campana  et  al.  2006;  Waxman  et  al.  2007 ;  Soderberg 
et  al.  2008),  although  alternative  scenarios  were  also  proposed 
that  invoked  the  presence  of  relativistic  jets  (Ghisellini  et  al. 
2007;  Mazzali  et  al.  2008).  The  UV/optical  counterparts  de¬ 
tected  in  co-aligned  Swift  UVOT  observations  were  well  fitted 
by  analytical  models  for  emission  from  the  expanding,  cooling 
envelope  of  a  progenitor  with  the  radius  of  a  W-R  star  (Campana 
et  al.  2006;  Soderberg  et  al.  2008;  Rabinak  &  Waxman  2010). 

The  wide-field  UV  survey  capabilities  of  the  Galaxy 
Evolution  Explorer  ( GALEX)  satellite  serendipitously  detected 
the  onset  within  7  hr  of  shock  breakout  of  a  UV-bright  plateau 
lasting  2  days  before  fading  away  from  Type  IIP  SN  SNLS- 
04D2dc  at  z  =  0.185,  two  weeks  before  its  discovery  in  the 
optical  Supernova  Legacy  Survey  (SNLS;  Gezari  et  al.  2008; 
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Figure  1.  Strong  Hcf  and  H/l  P-Cygni  features  in  SN  2010aq,  with  the  velocity  of  the  H/>  absorption  minimum  (which  shifts  from  14,000  km  s-1  to  6000  km  s_1 
between  0  =  10.7  and  53.2  days)  marked  with  a  dotted  line.  Also  shown  is  the  0  =  12.0  day  and  0  =  57.9  day  spectra  for  SN  2006bp  from  Dessart  et  al.  (2008) 
and  Quimby  et  al.  (2007)  which  have  been  multiplied  by  a  factor  of  0.013  and  0.002,  respectively,  and  have  an  1  \f>  absorption  minimum  at  10,160  km  s-1  and 
7,000  km  s-1,  respectively.  The  spectra  have  been  binned  in  wavelength  when  necessary  to  improve  the  S/N  and  offset  in  flux  for  clarity. 


Schawinski  et  al.  2008).  While  the  data  were  modeled  both 
numerically  and  analytically  with  thermal  emission  following 
shock  breakout  in  an  RSG  progenitor  with  Rt  ranging  from 
800  to  1 000  R@  and  a  variety  of  stellar  envelope  density  struc¬ 
tures  (Gezari  et  al.  2008;  Schawinski  et  al.  2008;  Tominaga  et  al. 
2009;  Rabinak  &  Waxman  2010),  the  low  signal-to-noise  ratio 
of  the  detections  (S/N  <  3)  and  the  lack  of  simultaneous  opti¬ 
cal  data  preclude  one  from  fitting  for  the  parameters  of  the  RSG 
progenitor  star. 

In  this  Letter,  we  present  SN  IIP  2010aq  at  z  =  0.0862 
discovered  in  joint  monitoring  observations  in  the  UV  by 
GALEX  Time  Domain  Survey  (TDS)  and  in  the  optical  by  the 
Pan-STARRSl  Medium  Deep  Survey  (PS1  MDS).  The  UV  and 
optical  observations  discover  the  SN  less  than  20.5  hr  after 
shock  breakout  in  the  SN  rest  frame,  and  the  high  S/N  of  the 
detections  enable  us  to  use  model  fits  to  the  early  light  curve  to 
place  meaningful  constraints  on  the  radius  of  the  RSG  progenitor 
star  for  the  first  time. 

2.  OBSERVATIONS 

From  UT  2010  February  4  to  March  22  (UT  dates  are  used 
throughout  this  Letter),  the  GALEX  TDS  monitored  the  7  deg2 
PS1  MDS  field  MD04  in  the  NUV  with  7  pointings,  each  1 
deg2,  with  a  cadence  of  2  days  and  with  a  limiting  magnitude 
of  wab  ~  23.  Gaps  in  the  regular  monitoring  occurred  when 
the  NUV  detector  was  shut  down  for  engineering  tests.  The 
GALEX  science  images  were  run  through  the  difference  imaging 
pipeline  photpipe  (Rest  et  al.  2005),  and  difference  image 
detections  were  filtered  through  a  series  of  selection  criteria  to 
remove  false  detections  due  to  image-subtraction  artifacts.  The 
remaining  detections  were  associated  between  epochs,  and  then 
cross-matched  with  the  PS1  SN  alerts  from  the  MDS  with  a 
matching  radius  of  5". 

GALEX  TDS  detected  a  variable  object  on  2010  February  12 
coincident  with  a  PS1  SN  candidate  discovered  at  R.A.= 
10h02m09‘‘:741,  decl.=  +01°14,00'.,94  (J2000)  in  daily  imaging 


cycling  through  the  g,  r,  i,  and  z  filters  starting  on  2010 
February  12.  The  SN  was  confirmed  with  spectroscopy  with 
the  4.2  m  William  Herschel  Telescope  (WHT)  with  ISIS  on 
2010  February  23.08  to  be  a  Type  IIP  SN  at  z  =  0.0862  from 
its  strong  Balmer  P-Cygni  features  (SN  2010aq;  Valenti  et  al. 
2010).  We  obtained  additional  epochs  of  spectroscopy  with 
the  Nordic  Optical  Telescope  (NOT)  with  ALFOSC  on  2010 
February  24.98,  the  6.5  m  Magellan  Telescope  with  IMACS  on 
2010  February  25.17,  and  with  the  6.5  m  MMT  with  Hectospec 
on  2010  April  10.20.  In  Figure  1,  we  show  the  Ha  and  H/3 
features  in  our  spectra,  and  the  velocity  drift  of  the  absorption 
minimum  of  H/3  with  time  from  the  recession  of  the  photosphere 
in  mass. 

In  order  to  avoid  systematic  errors  from  residuals  in  the 
difference  images,  the  GALEX  photometry  is  measured  on  the 
science  images  with  a  6"  radius  aperture  centered  on  the  host 
galaxy,  and  the  flux  of  the  SN  is  determined  by  subtracting 
the  flux  of  the  host  galaxy  measured  from  pre-SN  epochs 
(NUV  =  22.0  ±0.1  mag).  The  I  cs  error  is  measured  empirically 
from  the  dispersion  in  magnitudes  measured  for  sources  detected 
in  all  epochs.  The  offset  between  the  SN  and  the  host  galaxy 
centroid  is  1"5  (1  pixel),  for  which  we  estimate  an  upper  limit 
to  systematic  errors  due  to  centering  of  the  aperture  <0.1  mag. 
The  GALEX  photometry  is  in  the  AB  system,  has  been  aperture 
corrected  (Am  =  —0.23;  Morrissey  etal.  2007),  and  is  corrected 
for  Galactic  extinction  with  E(B  —  V)  =  0.024  mag. 

The  PS1  MDS  images  were  reduced  with  the  PS1  Image 
Processing  Pipeline  (Magnier  2006)  including  detrending,  cali¬ 
bration,  warping,  and  stacking.  Difference  imaging  and  tran¬ 
sient  alerting  were  performed  using  the  photpipe  pipeline 
on  the  stacked  PS1  images  (Rest  et  al.  2005).  The  SN  pho¬ 
tometry  is  measured  by  photpipe  using  forced  point-spread 
function  (PSF)  photometry  on  the  difference  images  at  the  po¬ 
sition  of  the  SN,  and  the  zeropoint  is  measured  from  com¬ 
parison  with  field  stars  in  the  SDSS  catalog.  The  PS1  pho¬ 
tometry  is  in  the  AB  system  and  is  corrected  for  Galactic 
extinction. 
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Table  1 

GALEX  and  PS  1  Photometry 


Filter 

Date  (MJD) 

4>  (days)3 

magb 

1(T 

NUV 

55235.42 

-2.8 

(22.69) 

NUV 

55237.54 

-0.8 

(22.68) 

NUV 

55239.32 

0.8 

21.06 

0.18 

NUV 

55241.52 

2.8 

20.64 

0.09 

NUV 

55249.66 

10.3 

21.69 

0.17 

NUV 

55251.44 

12.0 

21.68 

0.18 

NUV 

55259.52 

19.4 

(22.48) 

NUV 

55261.64 

21.4 

(22.37) 

NUV 

55266.98 

26.3 

(22.14) 

NUV 

55271.98 

30.9 

(22.55) 

NUV 

55277.39 

35.9 

(22.36) 

g 

55236.48 

-1.8 

(24.00) 

g 

55239.48 

1.0 

21.63 

0.08 

g 

55242.43 

3.7 

20.90 

0.03 

g 

55248.43 

9.2 

21.00 

0.03 

g 

55275.39 

34.0 

21.80 

0.07 

g 

55290.38 

47.8 

22.00 

0.24 

g 

55305.27 

61.5 

23.03 

0.32 

r 

55236.49 

-1.79 

(23.90) 

r 

55239.50 

1.0 

21.74 

0.07 

r 

55242.44 

3.7 

20.95 

0.02 

r 

55248.44 

9.2 

20.93 

0.03 

r 

55275.41 

34.0 

21.29 

0.05 

r 

55305.29 

61.6 

21.92 

0.11 

i 

55237.49 

-0.9 

(23.99) 

i 

55240.47 

1.9 

21.60 

0.05 

i 

55243.33 

4.5 

20.90 

0.03 

i 

55246.32 

7.3 

20.86 

0.04 

i 

55249.48 

10.2 

20.76 

0.03 

i 

55267.41 

26.7 

20.77 

0.05 

i 

55288.39 

46.0 

20.88 

0.05 

i 

55306.25 

62.4 

21.64 

0.09 

i 

55321.27 

76.3 

22.42 

0.14 

z 

55235.45 

-2.7 

(22.83) 

z 

55238.47 

0.0 

22.60 

0.22 

z 

55241.39 

2.7 

21.42 

0.07 

z 

55250.49 

11.1 

20.89 

0.06 

z 

55259.33 

19.2 

20.92 

0.05 

z 

55289.29 

46.8 

21.40 

0.07 

z 

55298.34 

55.2 

21.53 

0.10 

z 

55301.37 

57.9 

20.98 

0.16 

z 

55304.26 

60.6 

21.70 

0.08 

Notes. 

a  Days  since  MJD  55238.4  in  the  SN  rest  frame. 

b  AB  magnitudes  corrected  for  a  Galactic  extinction  of  E(B  —  V)  =  0.024  and 
the  host  galaxy  contribution  removed.  Upper  limits  are  given  in  parentheses. 


Table  1  lists  the  GALEX  and  PS  1  SN  photometry,  and  Figure  2 
shows  the  UV  and  optical  SN  light  curve.  We  define  the  phase 
of  the  SN  (0)  to  be  days  in  the  SN  rest  frame  since  the  time  of 
shock  breakout  (fs bo)-  Our  observed  constraint  on  ?sbo  is  the 
time  between  the  last  non-detection  with  GALEX  in  the  NUV 
on  MJD  55,237.54  and  the  first  detection  with  PS  1  in  the  z  band 
on  MJD  55,238.47,  which  is  a  window  of  20.5  hr  in  the  SN  rest 
frame.  We  further  constrain  fsBO  using  fits  of  shock  breakout 
models  to  the  early  UV /optical  light  curve  in  Section  4  to  be 
rSBO  =  MJD  55,238.42  ±  0.02. 


3.  HOST  GALAXY 

We  measure  the  host  galaxy  redshift  of  z  =  0.0862  from 
nebular  emission  lines  that  are  superimposed  on  the  SN  spec¬ 
trum  and  adopt  a  distance  modulus,  DM  =  38.0.  The  narrow 


Figure  2.  GALEX  NUV  and  PS1  g,  r,  i,  and  z  light  curve  of  SN  2010aq.  The 
flux  of  the  host  galaxy  has  been  removed,  and  the  magnitudes  are  in  the  AB 
system,  corrected  for  Galactic  extinction,  and  subtracted  by  the  distance  modulus 
(DM  =  38.0).  Solid  and  dashed  lines  show  the  light  curves  of  well-studied  Type 
IIP  SNe  2005cs  and  2006bp  from  Dessart  et  al.  (2008)  and  ^T-corrected  to  the 
observed  frame  of  SN  2010aq.  SN  2005cs  was  reddened  by  0.02  mag  and  SN 
2006bp  was  dereddened  by  0.24  mag  to  remain  with  a  net  internal  reddening  of 
E(B  —  V)  =  0.1 1  d=  0.05  mag  and  offset  by  —1.9  and  —0.3  mag,  respectively. 


Balmer  lines  have  a  Galactic-extinction  corrected  Hcr/H/i  ra¬ 
tio  of  3.28  ±  0.15,  corresponding  to  an  internal  extinction 
of  E(B  —  V)  =  0.12  ±  0.04  mag  (Calzetti  2001).  We  use 
the  semiempirical  relation  between  metallicity  and  the  mea¬ 
sured  nebular  emission  line  ratios,  [N  ii]/.6583/H«  <  0.1, 
[Oiii]X5007/H/3  =  3.6  ±  0.2,  and  log(£23  =  ([OnjA.3727  + 
[O  iii]X5007,  4959) /H/S)  =  0.89  ±  0.05,  to  measure  a  low  oxy¬ 
gen  abundance  of  12+log(0/H)  =  8.2-8.5  (McGaugh  1991; 
van  Zee  et  al.  1998;  Pettini  &  Pagel  2004),  which  is 
consistent  with  the  low  metallicity  expected  from  the 
luminosity-metallicity  relation  from  Tremonti  et  al.  (2004) 
given  the  dwarf  host  galaxy  of  SN  2010aq  with  an  internal 
extinction  and  /^corrected  rest-frame  absolute  magnitude  of 
Mb  =  —17.3  mag. 


4.  COMPARISON  WITH  ANALYTICAL  MODELS 

We  compare  the  early  (<f>  <  3  days)  UV /optical  light  curve 
to  analytical  models  for  shock  breakout  and  its  immediate 
aftermath  in  an  RSG  star,  the  progenitors  of  Type  IIP  SNe. 
The  analytical  model  presented  by  Nakar  &  Sari  (2010,  NS  10) 
follows  the  photon-gas  coupling  to  determine  the  location  of 
the  photosphere  and  takes  into  account  photon  diffusion  from 
inner  layers  below  the  photosphere.  The  analytical  relations 
presented  by  Rabinak  &  Waxman  (2010,  RW10)  neglect  the 
effects  of  photon  diffusion,  but  do  include  opacity  from  electron 
scattering  and  its  effect  on  the  emitted  spectrum.  Although 
the  pre-SN  stellar  structure  and  detailed  physics  are  better 
treated  in  numerical  radiation  hydrodynamics  models,  these 
models  are  computing-time  intensive,  and  analytical  models 
are  more  practical  to  use  for  fitting  for  the  parameters  of  the 
progenitor  star.  We  rely  on  the  analytical  scaling  relations  from 
NS10  to  constrain  R5 oo  and  E^/ where  £500  is  £*/500£q, 
M is  =  M/15M0,  and  £51  =  E / 1 05 1  ergs.  Here,  £  is  the  energy 
of  the  explosion,  and  M  refers  to  the  mass  of  the  ejecta  and  does 
not  include  the  1 .44  M0  from  the  neutron  star  core  or  mass  lost 
from  winds  during  the  lifetime  of  the  star. 

In  the  spherical  expansion  phase  (when  £ph  ~  i;ph?),  the 
NUV  light  curve  rises  to  a  maximum  when  the  ejecta  expands 
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Figure  3.  Zoom  in  of  the  early  UV/optical  light  curve  of  SN  2010aq  in 
comparison  to  the  best-fitting  analytical  model  from  Nakar  &  Sari  (2010, 
NS  10)  with  E51  =  1,  Mi 5  =  1,  and  R500  =  1  (thick  solid  line)  and  the 
analytical  model  from  Rabinak  &  Waxman  (2010,  RW10)  with  £51  =  1.44, 
M15  =  0.59,  and  £500  =  1-73  (dash-dotted  line).  The  models  have  been 
^-corrected  into  the  observed  frame  of  SN  2010aq  and  applied  with  a  reddening 
of  E(B  —  V)  =  0.03  mag.  The  NS  10  model  has  been  offset  by  — 1.5  mag. 


and  cools  until  the  peak  energy  of  the  emission  is  3 kT  — 
/(VnuvG  +  z)  or  when  T  —  1.94  eV.  Using  T(cf>)  for  an  RSG 
fromNSlO,  this  translates  to</>peak  =  2.2  days 
For  a  blue  supergiant  star  (BSG)  with  £  =  20-70 £0,  <j&peak  in 
the  NUV  occurs  much  earlier  (<j>peak  <  0.6  days),  which  is 
incompatible  with  the  observed  light  curve  of  SN  2010aq.  If 
we  use  our  observed  limits  on  fsBO  =  MJD  55,238.0  ±  0.5  and 
conservatively  estimate  <j>peak  to  be  within  1  day  of  the  second 
GALEX  detection  on  MJD  55,241.5,  then  in  the  SN  rest  frame 
(/>peak  =  3±1  d  and  this  results  in  M^0'23/?^8^^20  =  1.4±0.5. 

The  first  detection  in  the  NUV  ((pi)  on  MJD  55,239.30,  as 
well  as  on  MJD  55,239.46  in  the  g  and  r  bands,  places  a  tight 
constraint  on  T((p\ ).  In  the  first  couple  of  days  after  shock  break¬ 
out,  the  spectrum  of  the  emission  is  well  approximated  by  a 
diluted  blackbody  (Gezari  et  al.  2008;  Tominaga  et  al.  2009). 
Assuming  that  the  g  and  r  data  points  are  ~0.2  mag  fainter 
on  MJD  55,239.30,  the  UV/optical  photometry  is  fitted  with 
E(B  —  V)  —  0.11  ±  0.01  mag  and  7)*  =  2.7  ±  0.5  eV,  in 
excellent  agreement  with  E ( B  —  V)  derived  from  the  Balmer 
decrement.  Note  that  the  observed  diluted  blackbody  tempera¬ 
ture  is  hotter  than  the  photospheric  temperature  by  a  factor  of 
^1.2  as  a  result  of  electron  scattering  (RW10). 

In  Figure  3,  we  plot  the  NS10  model  for  an  RSG  with 
M15  =  R500  =  £51  =  1,  fitted  for  E(B  —  V)  —  0.03  mag  and 
rSB0  =MJD  55,238.42  ±0.02,  with  an  offset  of  —  1.5  mag  and  a 
/(-correction  to  transform  the  NS10  models  from  the  rest  frame 
to  the  observed  frame  of  SN  2010aq.  The  NS  10  model  is  in  good 
agreement  with  the  data  until  <p  =  3  days.  The  tight  constraint 
on  fsBO  is  from  the  fit  to  the  earliest  detection  in  the  z  band  on 
</>  =  0.03  ±  0.02  days  and  does  not  include  the  much  larger 
uncertainties  at  very  early  phases  (<p  <  1.5  days  for  an  RSG) 
introduced  by  the  simplistic  approximations  in  the  analytical 
models  (RW10).  While  the  NS  10  M15  =  £500  =  £51  =  1 
model  is  fitted  to  the  data  with  a  small  amount  of  reddening, 
if  we  require  that  E(B  —  V)  —  0.12  ±  0.04  mag,  then  one 
gets  the  same  fit  at  (pi  to  a  model  with  a  factor  of  1  .4/q2  times 
higher  temperature,  and  thus  with  'W|50J3£“8£°|I  1  =  1.4/°02. 
Given  the  offset  in  luminosity  required  to  fit  the  model  and 


L(<p)  for  an  RSG  from  NS  10,  one  gets  an  additional  constraint 
ofM^Rsoofi?!96^. 

When  combining  all  of  the  constraints  above,  an  extinction 
of  E(B  —  V)  <  0.12  mag  is  required  for  the  allowed  parameters 
from  (ppeak  to  converge  with  those  from  T(<p\).  If  we  assume 
the  best-fit  value  for  the  reddening  from  T(cp\)  of  E(B  —  V)  = 
0.11  ±  0.01  mag,  which  is  in  agreement  with  this  upper  limit, 
then  for  different  values  of  £51,  one  then  gets  the  following 
solutions:  (£51,  M15,  £50 0)  =  (0.5,  0.09-0.15,  1.0-1 .5),  (1.0, 
0.23-0.38, 1. 1-1.8),  and  (1.5, 0.37-0.62, 1.1-1. 8),  thus  placing 
the  strongest  constraint  on  the  radius  of  the  RSG  progenitor, 
R500  =  1.4  ±  0.4.  For  this  value  of  £*,  stellar  evolutionary 
models  of  core-collapse  SN  progenitor  stars  with  solar  to  half¬ 
solar  metallicities  have  M\$  ~  0.45-0.90  (Eldridge  &  Tout 
2004),  and  hydrodynamical  simulations  of  SNe  IIP  have  even 
larger  model  ejecta  masses  (M15  >  1;  Utrobin  &  Chugai  2008), 
thus  favoring  the  solutions  with  £51  =1.5. 

In  Figure  3,  we  also  show  a  /(-corrected  analytical  model 
with  a  reddening  of  E(B  —  V)  =  0.03  mag  from  RW10  with 
parameters  within  the  allowed  solutions  (£51  =  1.44,  M 15  = 
0.59,  and  £500  =  1.73)  that  has  a  luminosity  and  UV/optical 
colors  that  are  in  good  agreement  with  the  data  at  (pi,  but 
with  a  slower  temperature  evolution  that  is  in  slightly  worse 
agreement  with  the  data  at  (p  =  3  days.  If  we  again  assume 
E(B  —  V)  —  0.11  mag,  then  the  model  luminosity  increases  by 
a  factor  of  1.8,  resulting  in  a  similar  increase  in  £51/4/15  for 
values  of  £500  <  1-8. 

5.  COMPARISON  WITH  TYPE  IIP  SNe  2005cs 
AND  2006bp 

We  compare  the  late-time  (</>  >  3  d)  UV /optical  light  curve 
of  SN  2010aq  with  two  nearby  Type  IIP  SNe  well  studied  in 
the  UV  and  optical  with  the  Swift  satellite  and  ground-based 
photometry  and  spectroscopy,  SN  2005cs  (DM  =  29.7)  and 
SN  2006bp  (DM  =  31.2;  see  summary  in  Dessart  et  al.  2008). 
In  order  to  compare  the  light  curves  to  SN  2010aq,  we  use 
SN  IIP  template  spectra  to  measure  the  /(-correction  required 
to  transform  the  observations  into  the  GALEX  and  PS1  filters 
in  the  observed  frame  of  SN  2010aq  and  apply  a  net  internal 
reddening  of  £(£  —  V)  =  0.1 1  mag,  where  we  use  the  Cardelli 
et  al.  (1988)  law  and  the  internal  reddening  estimated  for  SN 
2005cs  (2006bp)  from  synthetic  spectral  fits  to  optical  spectra 
and  UV  photometry  of  E(B  —  V)  =  0.04  (0.40)  mag,  with  an 
accuracy  of  ±0.05  mag  from  Dessart  et  al.  (2008).  We  set  (p  =  0 
for  SN  2005cs  (2006bp)  to  the  jr-intercept  of  the  expanding 
photosphere  from  Dessart  et  al.  (2008)  of  MJD  53,547.3  (MJD 
53,832.9),  and  after  subtracting  the  distance  modulus,  shift  the 
light  curves  by  —1.9  mag  (—0.3  mag).  The  optical  light  curves 
are  in  good  agreement  between  the  SNe  in  Figure  2,  and  the 
similar  evolution  of  the  photospheric  velocity  with  time  of  SNe 
2010aq  and  2006bp  (see  comparison  in  Figure  1)  suggests  they 
have  a  similar  E/M,  since  i>ph  oc  +JE /M. 

Figure  2  shows  a  large  NUV  excess  of  SN  2010aq  compared 
to  SNe  IIP  2005cs  and  2006bp  at  (p  =  1 1-13  days  .  Given  that 
the  NUV  excess  is  more  prominent  at  (p  —  11-13  days  than  at 
(p  —  3  days  and  that  the  (p  =  10.7  day  spectrum  does  not  show 
signs  of  circumstellar  interaction,  it  is  most  likely  that  the  NUV 
excess  is  associated  with  the  effects  of  line  blanketing  by  metal 
ions  shortward  of  3000  A,  which  increases  with  time  when  the 
photospheric  temperature  of  the  SN  is  cool  enough  for  Fe  11  to 
be  present  (Dessart  &  Hillier  2005).  Indeed,  if  the  metallicity  of 
the  solar-metallicity  non-LTE  atmosphere-code  synthetic  spec¬ 
tra  fitted  to  SNe  2005cs  and  2006bp  from  Dessart  et  al.  (2008) 
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was  decreased,  there  would  be  a  significant  increase  in  the  emit¬ 
ted  NUV  flux  (see  Figure  9,  Dessart  &  Hillier  2005).  However, 
one  must  be  cautious  with  this  interpretation,  since  the  UV 
spectral  range  is  difficult  to  model.  Even  in  the  case  of  SN 
1987A,  where  the  metallicity,  BSG  progenitor,  and  explosion 
energy  of  the  SN  are  well  constrained,  the  most  sophisticated 
time-dependent  non-LTE  radiative  transfer  models,  while  in  ex¬ 
cellent  agreement  with  the  optical  observations  (at  the  level  of 
10%),  are  in  much  worse  agreement  in  the  UV  (only  to  the 
30%-50%  level;  Dessart  &  Hillier  2010).  However,  given  the 
low  metallicity  of  the  dwarf  host  galaxy  of  SN  2010aq,  a  sub¬ 
solar  metallicity  progenitor  star  is  the  most  natural  explanation 
for  the  UV  excess  at  late  times. 

6.  CONCLUSIONS 

Although  the  duration  of  the  radiative  precursor  from  shock 
breakout  is  too  short  to  be  resolved  by  rolling  SN  surveys,  the 
thermal  expansion  phase  in  the  days  following  can  be  studied 
in  detail.  In  this  Letter,  we  have  shown  that  with  daily  cadence 
UV  and  optical  imaging  we  can  use  analytical  models  to  fit 
for  the  radius  of  the  pre-SN  progenitor  star,  and  even  use  the 
evolution  of  the  UV  light  curve  to  directly  probe  its  metallicity. 
UV  observations  are  also  critical  for  constraining  the  time  of 
shock  breakout  in  distant  SNe  (fpeak-NBO  is  only  2  days  in  the 
NUV  compared  to  2  weeks  in  the  r  band)  and  can  improve 
the  accuracy  of  using  Type  IIP  SNe  as  distance  indicators 
(Schmidt  et  al.  1992;  Dessart  et  al.  2008;  Kasen  &  Woosley 
2009). 

Based  on  the  detection  efficiency  of  GALEX  for  detecting 
UV  emission  from  Type  II  SNe  measured  from  serendipitous 
observations  of  SNLS  SN  candidates,  Gezari  et  al.  (2008) 
predict  a  detection  rate  for  GALEX  TDS,  which  monitors  one 
PS1  MD  field  per  month,  of  ^  1  month-1.  Given  the  actual 
detections  of  two  SNe  IIP  in  effectively  2  months  of  operations 
of  GALEX  TDS  in  coordination  with  PS  1 ,  luminous  SN  2009kf 
at  z  =  0.182  (Botticella  et  al.  2010)  and  SN  2010aq,  our 
detection  rate  appears  to  be  right  on  the  mark.  Continued  GALEX 
TDS  and  PS1  MDS  monitoring  will  yield  an  exciting  data  set 
to  probe  the  diversity  of  progenitors  of  core-collapse  SNe  and 
approach  a  better  understanding  of  the  mapping  from  stellar 
evolution  to  stellar  death. 
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